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We derive the exact equation of motion for a vortex in two- and three-dimensional non-relativistic 
systems governed by the Ginzburg-Landau equation with complex coefficients. The velocity is given 
in terms of local gradients of the magnitude and phase of the complex field and is exact also for 
arbitrarily small inter-vortex distances. The results for vortices in a superfluid or a superconductor 
are recovered. 
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Vortices are found in a variety of physical systems. Ac- 
cordingly, the study of these intriguing collective excita- 
tions attracts widespread attention among the physics 
community. Examples of vortices often studied are hy- 
drodynamic vortices, vortices in superfluids, in supercon- 
ductors and in nematic crystals, and cosmic strings ||],||. 
An important goal is to clarify the mechanisms by which 
vortices are created and the details of their motion sub- 
ject to local interactions, such as crossing, merging and 
intercommutation, as well as long-range forces. These 
issues have recently been addressed in the context of rel- 
ativistic scalar field theories HQ- 

In this Letter we present an analytic derivation of the 
exact equation of motion for a vortex in a non-relativistic 
dissipative system. The system we study is one modelled 
by the extensively studied Ginzburg-Landau equa- 
tion with complex coefficients (CGL) 

£ t A = P(A,A*)A + b\7 2 A (1) 

where A = \A\ exp(iS) is a complex field, the function P 
is given by P(A, A*) = /.i - a\A\ 2 , and a,b,fi£ C. By a 
suitable rescaling of time, space, and A, the number of 
(real) adjustable parameters in the coefficients of eq. (Q) 
may be brought down to two, as is often done. However, 
we shall keep eq. (Q) unsealed for clarity. We study the 
equation in two and three spatial dimensions. 

The reason for selecting the CGL equation is two-fold. 
Firstly, it is a relatively simple partial differential equa- 
tion; yet it exhibits the principal features of more compli- 
cated oscillatory systems. A prime example of such sys- 
tems are reaction-diffusion systems, such as the chemical 
oscillatory Belousov-Zhabotinsky reaction |Hj,|ll|] • 

Secondly, the CGL equation contains a number of in- 
teresting special cases. When a, b, and [i are purely imag- 
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FIG. 1. One-armed spiral vortex of the CGL field 
A = |yl|e' s in a two-dimensional system. The height of 
the surface depicts the magnitude \A\. The spiraling curves 
are contour lines of the phase S (isophase lines). The phase 
change between two thin lines is n/2. 

inary the CGL equation coincides with the non-linear 
Schrodinger equation. The latter equation describes the 
quantum dynamics of superfluid 4 He and is known in 
that context as the Ginzburg-Pitaevskii-Gross equation 
(GPG) fl^| . Furthermore, by employing the Madelung 
transformation ]T^ | the non-linear Schrodinger equation 
also transforms into the hydrodynamic equations for an 
inviscid fluid (the Euler equations). In both cases \A\ 2 
corresponds to the (super)fluid mass density and VS* is 
proportional to the velocity of the (super)fluid. We stress 
that the GPG case is special because it describes a conser- 
vative system and the vortex motion can be derived from 
a Lagrangian. The CGL equation, on the other hand, de- 
scribes a dissipative system and one is compelled to pur- 
sue a direct derivation of the vortex equation of motion, 
as we do here. 

Equation ([!]) permits solutions in which A has phase 
singularities (defects). In two space dimensions these 
are isolated points around which the phase S changes 
by multiples of 27r. At the same points the magnitude 
\A\ vanishes, so that the complex field A remains single 
valued, see figure In the vicinity of a defect the phase 
is of the form S = ip( r ) + ^(v 3 — tjJ ^) i n polar coordinates 
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(r, if) |L4j . For a constant phase S this is the equation 
for n|-armed spirals rotating at an angular frequency iv. 
In three dimensions the defects become one-dimensional 
strings, or filaments, and the spirals generalize to scroll 
waves 0,^5) which look like sheets wound around a fila- 
ment. The filaments may be closed or open (in which case 
they end on the system boundary) and of arbitrary shape. 
We shall call a solution with one defect or filament (in 
two or three dimensions) a spiral vortex, in analogy with 
the (non-spiral, i/j'(r) — 0) vortex solution of the GPG 
equation, which describes the circulation of the super- 
fluid around strings of normal-phase fluid. The integer n 
is the winding number of the vortex and is a topologically 
conserved quantity in two dimensions but not in three. 
The core of the vortex is the region where the magnitude 
| .A | deviates significantly from its asymptotic value, see 
figure 0. 

The evolution of a system with (spiral) vortices may 
be described in terms of the motion of the defects, or 
filaments, along with values of the fields \A\ and S at 
positions away from the defects or filaments. Such a sep- 
aration into collective coordinates and field variables is 
non-trivial, and the present work comprises the first ex- 
act treatment of this kind for a dissipative system. The 
motion of a vortex is affected by modifications in the field 
A due to the presence of other vortices or system bound- 
aries. If the vortices are assumed to form a dilute system, 
i.e. one where the defects are well separated, the influence 
of variations in the magnitude \A\ of the complex field 
may be neglected, since \A\ will assume its asymptotic 
value at distances much smaller than the inter-defect dis- 
tance |lq| . Under this assumption, the interaction be- 
tween vortices can be described entirely by the phase S. 
In this approximation Rica and Tirapegui || (and in a 
slightly different form also Ref . || ) have derived the equa- 
tion of motion in two space dimensions for the position 
of the kth defect X^(t) in terms of the portion of the 
phase S due to other defects, 9^ k \x) = S — rikPk: where 
tanc/?fc = (y — Yk)/(x — Xk)- Their result (for |nfe| = 1 
and b R — 1, but here generalized to any value of and 
b R ) is 



X, 



dXj 
dt 



26 7 VffW _ 2b R ^z x V0 W , (2) 
\n k \ 



where b R = Re b, b : = Im b, and z — x x y is normal to 
the plane. The first term, proportional to the gradient, is 
that found by Fetter in the GPG limit corresponding 
to b R = 0, 6j = h/ (2m)and states that the vortex moves 
with the local superfluid velocity. The second term is the 
perpendicular Peach-Koehler term jl8| first found in this 
context by Kawasaki [|l9| . 

When the system of spiral vortices cannot be approx- 
imated by a dilute system the expression (||) for the de- 
fect velocity is no longer valid but will acquire additional 
terms. We shall take a completely general approach in 
which the amplitude \A\ is allowed to vary. This will en- 
able us to determine the exact motion of a defect also 



when another defect is located an arbitrarily small dis- 
tance away, i.e. even when the vortex cores overlap. It 
will also provide the exact motion of a defect which is 
arbitrarily near a system boundary. For filaments in a 
three-dimensional system our treatment will furthermore 
correctly incorporate interactions with other segments of 
the same filament. 

The corresponding problem for a relativistic scalar field 
theory was solved by Ben-Ya'acov Q . His derivation was 
based strictly on a covariant world-sheet formalism that 
cannot be applied to a non-relativistic theory. For the 
CGL equation one must therefore resort to other meth- 
ods. 

Let us consider the general motion of vortices in three 
space dimensions. The motion in a two-dimensional sys- 
tem can be found from the three-dimensional problem as 
the special case of straight, aligned vortices. 

We may generalize eq. ([!]) by admitting any continu- 
ous function P(A,A*). The details of P do not enter 
the derivation. The field A is zero on a collection of 
one-dimensional strings which are the filaments. Let the 
position of the filament T of a vortex be given at time t by 
X(s,t), where s is the arclength coordinate along T. We 
define a local coordinate system along the string as fol- 
lows Q . At each point along the string the unit tangent 
vector T — dX/ds, the unit normal vector N, and the 
binormal vector B = T x TV form an orthonormal frame 
so that any position x in a neighborhood of the string can 
be expressed as x = X(s,t) + xN(s,t) + yB(s,t). The 
coordinate representation (s,x,y) is unique for x < 1/k 
but becomes singular when x reaches or exceeds the ra- 
dius of curvature 1/k. 

Along the string, the transport of the unit vectors is 
given by the Frenet-Serret equations pffl 
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where k is the curvature and r is the torsion of the string. 
Let us further introduce the local polar coordinates r, ip 
defined by x = r cos ip, y = r simp. In terms of these 
coordinates, the gradient and Laplacian take the forms 



V = TH + r^- 

or 

dr 2 

K 



H 2 
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cos ip 
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sin ip 
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where 



H _ 1 d r d 

1 — nr cos ip \ ds dip 



(4) 



(5) 



(6) 



We now proceed to find the velocity X(s,t) of the fila- 
ment r. Because this string of zeros of the function A 
has no transverse extension and is a feature of a solu- 
tion of an underlying local field theory, its motion should 
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be determined from the behavior of the fields |A| and 
S in an infinitesimal neighborhood of the filament. It 
will be sufficient to study the fields within a distance 
e -C min(<i, where d is the shortest distance to an- 
other string segment pl[ . This condition ensures unique- 
ness of the coordinate representation. 

The phase field S is multi-valued and satisfies 
S(s,r,tp + 2n\t) - S(s,r,tp;t) — n 2n for < r < e. 
Let us therefore split S — x + # m such a way that \ 
contains all multi-valued contributions to the phase and 
depends on time only through the position of the fila- 
ment r. For a straight (or two-dimensional) isolated vor- 
tex one may choose x — nL P- A consistent description 
of the multi- valued phase of an arbitrarily shaped vortex 
filament requires, however, a global realization such as 
the Biot-Savart integral, 



_ n f , _, x — X 



(7) 



This expression is known to contain logarithmic diver- 
gencies as r — > 0, as well as functions of the azimuthal 
angle ip that are multi- valued at r — We therefore 
absorb in x any part of S that is non-diffcrcntiable at 
r = 0. Similarly, we may write \A\ = Rw, where lni? 
depends on the filament position and contains all contri- 
butions to In \A\ that are non-differentiable at r = 0. For 
a straight isolated vortex one may choose R = r'"' . Thus 
9 and Inw are differentiable and it follows that the time 
derivatives 9 and w are finite for r < e. We remark that 
the choice of x an d R is not unique, since S and \A\ are 
invariant under two independent local symmetries 



X -> X + 
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R^Rf, w -> wf' 1 , (8) 



where 6 and In / are differentiable. 

With these definitions the real and imaginary parts of 
equation (|j) lead to the two equations 



d_ 

dt 



where 



(lni? 



dt 



Inw) = Re(P) 
> + X )= Im(P) 



b R Qi 
■hQi- 



b R Q2 



(9) 
(10) 



Qi = V 2 lni? + V 2 lnw 

+ (Vlni? + Vlnw) 2 - (Vx + V6>) 2 , 
Q2 = V 2 x + V 2 6» + 2(Vlni? + Vlnw) • (Vx + V6>). 

The time derivative <i/cft in eqs. @ and (|l0|), which is 
to be evaluated in the lab frame, is related to the time 
derivative d/dt in the moving reference frame of the local 
segment of the filament by d/dt = —X ■ V + d/dt. 

In order to include logarithmic divergencies as well as 
multi-valuedness as r —* 0, we write V% = /if + [n/r + 
f 2 )(p + AiT and Vlni? = {\n\/r + f 3 )r + U<p + A 2 T, 
where 

fi(r,ip, s,t) = gi(<p, s,t) + hi(<p, s,t) Inzer + O(r) (11) 



and 0(r) denotes any terms that vanish as r — > 0. It 
can be easily confirmed from this equation (as well as ar- 
gued on general grounds) that dx/dt, dx/ds, d(ln R)/dt, 
d(\nR)/ds, Ai and A2 have well-defined finite limits as 
r — > 0. We require that Vx and Vlni? be integrable, 
and that they satisfy the following condition near the 
filament: 



Vx 



-T x VlnR = C(s,t) + 0(r) 



(12) 



The arbitrary vector C corresponds to a choice of gauge 
in eq. (JsJ) . In the symmetric gauge R — r'"', x = ^V 3 for 
a straight (or two-dimensional) isolated vortex we have 
C = 0. 

Since 9 and Inw are differentiable, the singularities of 
Vx and Vlni? at r = must satisfy eqs. (||) and (|10| ) 
order by order. This last condition together with eq. (113) , 
leads to the coupled non-linear system 



V In R ■ u + b R qi — b z q2 = regular 
Vx • u + b,qi + b n q 2 = regular 



(13) 



where q x = V 2 lni? + (Vlni?) 2 - (Vx) 2 , q-i = V 2 x + 
2Vlni? • Vx and u = X + 2b R (Vlnw + AT x V6») - 
26,(V0 - AT x Vlntu). 

Cancellation of terms of order in eq. ( |l3|) leads 
to two equations for the perpendicular components of 
u. The integrability condition provides four first-order 
differential equations relating the functions gi and hi, 
and together with four algebraic relations resulting from 
eq. (|l^) the system can be solved in terms of four con- 
stants of integration. The perpendicular components of 
u arc then uniquely determined in terms of C. Further- 
more, the singular terms of order r~ 1 lnKr in eq. dl3| ) 
cancel. It is always possible to set the tangential veloc- 
ity, which is void of physical meaning, to zero by a time- 
dependent reparametrization s — > s(i). In the language 
of relativistic string theory, this is referred to as world- 
sheet reparametrization invariance. The exact result for 
the velocity of the vortex filament is 



X = bj\ k— B + 2(V ± 9 + C ± ) -2— T x Vlnw 
* n n 



+ b R ( kN - 2V ± Inw - 2^-T x (V0 + C) ) , (14) 
V W / 

where ( )j_ = T x [T x ( )] and the fields on the right- 
hand side are to be evaluated at the filament position 
X(s,t). The exact two-dimensional result is obtained as 

K — ► 0. 

The value of X is independent of the choice of gauge 
for i? and x- Indeed, substituting C from eq. ([[2]) into 
eq. (n4j) we obtain the manifestly invariant expression 



X = lim 

r^0 



77 77 

bA k—B + 2V ± S - 2—T x V In \A I 

n\ \n\ 

bJ kN - 2V ± In Ml - 2^-T x VS ' 
\n\ 1 



(15) 
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in which the filament velocity is written in terms of 
gradients of the magnitude and phase of the original 
complex field A. Let us define the complex velocity 
Z = (TV + iB) ■ X and express the derivatives in eq. (|15| ) 
in terms of z = x + iy and its conjugate z* . Then a quite 
beautiful result emerges: 

Z = b (-4#r hiA(z, z*) + k) , n > 1 

(16) 

Z* =b(-4§^lnA(z,z*) + K) , n < -1 , 

where the right-hand side is to be evaluated at z = z* = 0. 

The results are to be interpreted as follows: The ve- 
locity of the central filament of a vortex gets contri- 
butions from the curvature n of the filament and from 
local gradients of the magnitude \A\ and phase S of 
the complex field. A cylindrically symmetric solution 
A = p(r) exTp[i(ip(r)+n(ip— tot))], for which p = \A\ ~ r'™ 
and ip'(0) = 0, contributes nothing to the velocity and 
corresponds to a straight (or two-dimensional) isolated 
vortex at rest with respect to the lab frame. Non-zero 
gradient contributions appear as a result of deviations 
from cylindrical symmetry in \A\ and 5*. In a symmet- 
ric gauge with C = 0, these deviations are represented 
by w and 9. The asymmetries arise from the presence 
of other vortices, system boundaries, or (in three dimen- 
sions) other segments of the same filament, causing the 
vortex to move. 

In the C = gauge the expression (Q) reproduces a 
variety of results obtained previously for special cases. 
For k = and V In w « it reduces to eq. (0) corre- 
sponding to a two-dimensional dilute system [pL In the 
GPG limit b a = the expression ( |l4|) coincides with that 
derived by Lee [Q, who used a different method to find 
the velocity. For bj = 0, eq. ([!]) describes the non-linear 
diffusion of two fluid components with identical diffusion 
constants. In this limit the contribution to X from cur- 
vature, b R nN, agrees with the result of Ref. Jl5| . 

The expressions (|l^)-(^6|) for the velocity are exact 
also for an arbitrarily small distance between filaments. 
This makes the formulation well suited for theoretical or 
numerical investigations of local vortex interactions, such 
as crossing, merging and intercommutation, in which the 
vortex cores overlap We caution that the GPG 

equation does not provide a realistic model for the core 
of a superfluid vortex, since there the core width is com- 
parable to interatomic distances. For magnetic flux vor- 
tices in a superconductor, however, the core width is 
much larger and a classical description is justified. Such 
vortices are solutions of eq. (|l]) with the substitution 
V — > V + 2ieA/ (he), where A is the vector potential and 
2e is the charge of a Cooper pair. The corresponding 
filament velocity is easily obtained by adding 2eA/(hc) 
to V0 in eq. Q or to VS in eq. (jl|) g§. 

In summary, we have derived the exact equation of 
motion for a vortex in a large class of models of a 
non-relativistic complex field described by the complex 
Ginzburg-Landau equation (|lj) with an arbitrary, con- 
tinuous function P(A,A*). The velocity is expressed in 



terms of local gradients of the magnitude and phase of the 
complex field A. The result agrees with that of Ref. Q 
(our eq. ||) in the case of a dilute two-dimensional system 
of vortices, but for the general non-dilute case in two and 
three dimensions we find additional contributions to the 
velocity corresponding to the asymmetry of the magni- 
tude | A around the vortex. 
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